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Abstract. An asymptotic solution of the system of Schwinger-Dyson equations for 
four-dimensional Euclidean scalar field theory with interaction ^((j)*(j)) 2 is obtained. 
For A > A cr = 167T 2 the two-particle amplitude has the pathology-free asymptotic 
behavior at large momenta. For A < X cr the amplitude possesses Landau-type sin- 
gularity. 
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1 Introduction 

The problem of asymptotic behavior at large momenta (or, at short distances) is one 
of oldest problem of quantum field theory. In general, this problem has been solved 
only for the asymptotically-free models. A solution of the problem of asymptotic 
behavior for other strictly renormalized models - quantum electrodynamics, self- 
interacting scalar field, Yukawa interaction - requires to go out the framework of 
weak-coupling approximation. 

The first attempt to define the asymptotic behavior in quantum electrodynam- 
ics was made by Landau and coworkers in the 1950s [1]. The result was unrehearsed: 
the photon propagator included the non-physical pole in Euclidean region of mo- 
menta. Then similar poles were indicated in other strictly renormalized models: 
self-interacting scalar field and Yukawa interaction. The presence of such singular- 
ities in Euclidean region violates general principles of the local field theory and is 
a serious problem for these models. 

A method of removing these non-physical singularities was proposed soon in 
works [2] . This method is based on the application of Kallen-Lehmann representa- 
tion to restore the correct analytical structure and remains the basic tool to solve 
the problem of the non-physical singularities to our time. However, the absence of 
dynamical foundation enforced to search a more detailed investigation for the role 
of these singularities. 

Further development of the quantum field theory has demonstrated that these 
non-physical singularities arise practically inevitably in the framework of any known 
non-perturbative methods: at the renormalization- group summation, in the frame- 
works of l/A^-expansion and mean-field expansion, etc. 2 

A widespread opinion is formulated as a mathematical inconsistency of the 
quantum field models that are not asymptotically free. There is rigorous theorem 
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[4] that the four-dimensional scalar field theory with 4 interaction on the lattice 
does not have an interacting continuum theory as its limit for zero lattice spacing, 
i.e. the theory is trivial. Weinberg argues [5], however, that this argument is 
inconclusive due to an uncertainty of the continuous limit in this model. 

Today the situation with triviality of 4> 4 theory is vague as before. Recent papers 
(see, e.g., [6], [7], [8] ) in this topic maintain incompatible statements. While works 
[7] confirm the triviality scenario, works [6], [8] state non-trivial behavior in the 
strong-coupling limit. 

The asymptotic short-distance region in the models without asymptotic freedom 
is the region of strong (or, more exactly, non-weak) coupling, and it is the main 
difficulty in its investigation. It seems that the above-mentioned standard methods 
are too tethered to the weak-coupling region and are not enough meaningful to 
inform us about behavior at short distances in these models. In this paper the new 
approximation for the </> 4 theory is investigated. This approximation is based on 
a system of Schwinger-Dyson equations (SDEs) for the propagator and the two- 
particle function. The first equation is the exact SDE for the propagator, and the 
second equation is a truncated SDE for the two- particle function ("two-particle 
approximation"). This approximation is considered as a first non-trivial step of 
the sequence of general n-particle approximations, which tends to the exact infinite 
system of SDEs at n — > oo. An investigation of some truncation of SDEs is not 
a news, of course. A new step is the consideration of a system of these equations 
instead of the usual consideration of some single truncated equation. 

A structure of the paper is as follows: in section 2, the necessary notations 
and definitions are given; SDE for the generating functional of Schwinger functions 
is introduced in the formalism of a bilocal source. Using of the bilocal source is 
an essential point of the construction. We consider using this presumably as a 
convenient choice of the functional variable. In particular, this variable is very 
convenient for the construction of the mean-field expansion (MFE) by method of 
work [11], which is presented in section 2. The existence of the Landau pole in the 
two-particle amplitude of the leading approximation of MFE is also demonstrated 
in this section. 

In section 3, a general construction of the approximation scheme for the sys- 
tem of SDEs is given. In section 4, an alternative view to the approximation of 
preceding section is given. The system of equations for the propagator and the 
two-particle function is considered as a basic approximation for the construction 
of the modification of MFE of section 2. In section 5, the renormalization of the 
system of equations is made and some supplement simplifications are discussed. 

In section 6, the asymptotic solution of the equation for the two-particle am- 
plitude at large momenta is presented. In section 7, the asymptotic behavior of the 
amplitude at large momenta is discussed. The amplitude in this model possesses 
a self-consistent behavior (as a constant plus a decreasing oscillating term) at the 
values of the renormalized coupling A € (A cr , 2A cr ), where X cr = 16-7T 2 . At A < \ cr 
the amplitude has some Landau-type singularity in the pre-asymptotic region. At 
A > 2A cr the method of solution is, strictly speaking, cannot be applied due to the 
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negative value of the field-renormalization constant, but if one assumes formally 
such values, the results of the permitted interval can be continued to the region of 
strong coupling. Conclusions are presented in section 8. 



2 Preliminaries. SDE for the generating functional, 
mean-field expansion and Landau pole 

Consider the theory of a complex scalar field <j>(x) in a four-dimensional Euclidean 
space (x G E4) with the Lagrangian 

C = -d^*d^^ml ( p^-^*<f>) 2 (1) 

in the symmetric phase (toq > 0, A > 0). The generating functional of 2n-point 
(n-particle) Schwinger functions can be written as a functional integral 



G = y D((j), 0*)exp|y dxC(x) — j dxdycj)* (x)rj(x,'< 



y)4>(y)\, (2) 



where rj(x, y) is a bilocal source. The nth derivative of G over 77 with the source 
being switched off is the n-particle Schwinger function. The propagator of the field 
(f> is 

SG 

A(x-y) =< 4>{x)4>* (y) >- 



6r](y,x) 



r)=0 



the two-particle function is the second derivative of G, etc. 3 

The translational invariance of the functional integration measure leads to re- 
lation 

/ Wr) ^^ * (2;)exp {/ dxC ^- j dxdy<t>*(xMx,y)cl>(y)} = 0, 

which can be rewritten as the functional-differential SDE for generating functional 
G: 

{ml - dl) ^xj + / dy ^ y ^s£y7) + 5{x - y)G = \(JiL,*r (3) 

To construct the mean-field expansion for the generating functional we consider as 
a leading approximation the following equation 



3 A formalism of the bilocal source in the quantum field theory was first elaborated by 
Dahmen and Jona-Lasinio [9]. For further development and references see, e.g., [10] 
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whose solution is (in operator notations) 

G<°> - exp{-A • rj\, 
where A (p) = (m 2 + p 2 )^ 1 and m 2 = to 2 , + AA (x = 0). Here A (x = 0) = 

,4 

J A ^ 27r ^ 4 Ao(p), and label A means that some regularization implied. 

The leading approximation generates linear iteration scheme for SDE (3) 

G = G (0) + G (1) + • • • + G (n) + • • • , 

where 

A , , u ; : -(ml-dl)—- —6(x-y)G^ = / dy lV {x, yi )— -. (5) 

OT){x 1 x)dri{y 1 x) or]{y,x) J ot](y,yi) 

An idea of this iterative scheme is as follows: we shall consider "equation with 
constant coefficients" (4) as a leading approximation, i.e., equation (3) with the 
second term omitted. This term contains the source rj manifestly. The Schwinger 
functions are the derivatives of G(rj) in zero and only the behavior of G near r\ = 
is essential, therefore such an approximation seems to be acceptable. The term 
omitted contains the source and should be treated as a perturbation. 

The general solution of equation (5) is the functional G^ — P 2n G^\ where 
P2n is a polynomial of 2nth degree on source rj. Therefore at the nth step the 
computation of Schwinger functions reduces to solving a closed system of 2n linear 
integral equations. 

The first-step generating functional G^ reads 

g( 1 ) = Qg 2 V-a 1 -„)g(°). 

Here Gi is the leading-order two-particle function and Ai is the next-to-the-leading- 
order correction to the propagator. 

The iteration equation at n = 1 gives the equation for the two-particle function 

(m 2 -^X? 2 (* V y , ) =5(x-y')A (x>-y)-\A (x-y)G 2 ( * * ) , (6) 

whose solution is 

G 2 ( I V y , ^=A {x-y')A {x'-y)- 
- j dx 1 dx 2 A {x - xi)A (a;' - x 2 )f(x 1 - x 2 )A (x 1 - y)A (x 2 - y'), (7) 
where in the momentum space 



and 

Lo(p) = J A j^ Aoip + q)Aoiq) (9) 

is the single scalar loop. 

Note, that first term in formula (7) is the missed disconnected part of the two- 
particle function of a leading approximation. Hence, the connected structure of the 
two-particle function is restored at the first step of iterations. Such a peculiarity 
of the iteration scheme is originated by the bilocal source and is not something 
exceptional: as is well-known, the similar phenomenon appears also in constructing 
l/A^-expansion in the bilocal source formalism. The crossing properties of two- 
particle function in such iteration schemes are also restored stage-by-stage when 
the next terms in the expansion are taken into account (a discussion of this topic 
see in work [12]). 

All above formulae for the propagator A and the amplitude / contain divergent 
integrals and should be renormalized. In correspondence with the standard recipe 
we introduce the renormalized Lagrangian 

£ = -0 M 0*0 M 0-mW-^*0) 2 , (10) 

where <j>, m and A are the renormalized field, mass and coupling, and add the 
counter-terms 

AC = -(z - 1)8^*8^ - 5m 2 cf,*cf> - (z x - l)^(</>*^) 2 , (11) 

which absorb the divergences. 

Full Lagrangian £ b = C + AC can be written as 

C b = -d^* b d^ b - m 2 ^ b - y (<^ h ) 2 , (12) 



2 to 2 + Sm 2 x \z\ 



where 

4> b = y/z<j>, ml = ~ r """ . ^ b = (13) 
z z l 

Then all above calculations are reproduced with Lagrangian (12), i.e., with the 
replacement m 2 , — > m 2 , A — > X b , A — > A b , G2 — > G2b, f — > fb etc., and the nor- 
malization conditions are imposed on the renormalized propagator A and amplitude 
/. For the easement of the following calculations we choose the normalization point 
at zero momenta. 

The normalization conditions for the propagator A(p 2 ) = z~ 1 Ab(p 2 ) are 

A- 1 (0)=m 2 , (14) 



- 

dp 



2 A | P 2 =o = 1. (15) 
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These conditions define the mass-renormalization counter-term Sm 2 and the field- 
renormalization constant z. It is easy to see that z = 1 in the case, and, conse- 
quently, A -1 = A^ 1 = m 2 +p 2 , G2 — G2b- Thus, the amplitude / = is 

The normalization condition for the amplitude is 

/(0) = A. (17) 

This condition defines the coupling-renormalization constant z\ and the renormal- 
ized amplitude 

- YrkwY (18) 



where 



L r (p)=Lo(p)-Lo(0) = -^-5 / (fa!og(l + «(l-«)-4) (19) 

iD7r j m 

is the renormalized loop. 

As it follows from equations (18) and (19) the renormalized amplitude / pos- 
sesses a non-physical singularity (Landau pole) in the point p 2 = M 2 , where M\ is 
a solution of the equation 

1 + XL r (M 2 ) = 0. 

This equation has a solution at any positive A since function L r covers all values 
from to —00. At p 2 — > 00 

1 , P 2 
wtt z m z 

and Ml = mexp{^-}. As was yet noted in the introduction, the same Landau 
pole arises at the calculations of renormalized amplitude by other methods: in the 
frameworks of 1 /TV-expansion and renormalization-group summation. 



3 The system of SDEs and two-particle approxi- 
mation 



From SDE (3) one can go to the functional Z = log G. The SDE for Z is 

5 2 Z 



6rj(y,x) 



$v(y,yi) 



-+- 



5Z 



sz 



5ri(x,x)8r)(y,x) 6r}(x,x) drj(y,x) 
(20) 
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The system of SDEs is an infinite set of equations for n-particle functions Z n = 
5 n Z/5r] n \ ri =o- The first SDE is simply equation (20) with the source being switched 
off: 

(m 2 -%)A(x-y) = 6(x-y)-\Z 2 ( X x (21) 

where m 2 = nip + \A(x = 0) and A = — 5Z/5rj\ v=li is the propagator (or, the 
one-particle function). The second SDE is the derivative of (20) with the source 
being switched off: 



(m 2 - d 2 x )Z 2 



X 


y 


x' 


y' 




X 


< 






x 1 




(22) 



and so on. The nth SDE is the (n — l)th derivative of SDE (20) with the source 
being switched off and includes a set of functions from one-particle function A to 
(n + l)-particlc function Z n+ \. We call " the n-particle approximation of the system 
of SDEs" the system of n SDEs in which the first n — 1 equations are exact and 
nth SDE is truncated by omitting the (n + l)-particle function. It is evident that 
the sequence of such approximations goes to the exact set of SDEs at n — > oo. The 
one-particle approximation is simply equation (21) without Z 2 . This approximation 
has a trivial solution which is a free propagator. The two-particle approximation 
is a system of equation (21) and equation (22) without Z 3 : 



(m 2 -dl)Z 2 ^ X x , y^=S(x-y')A(x'-y)-XA(x-y)Z 2 ^ l 



(23) 



which includes A and two-particle function Z 2 . This non-linear system will be the 
object of present investigation. 

The idea of the present approximation scheme is very simple and natural. How- 
ever, the calculations became more and more complicated at each following stage: 
e.g., the three-particle approximation is a system of three non-linear equations for 
the propagator, the two-particle function and the three-particle function. In follow- 
ing section some other scheme will be considered, which leads to same equations 
(21) and (23), but is more realistic in the calculational sense. 



4 The particular solution of SDE for the generat- 
ing functional and modified mean-field expan- 
sion 

Another view to the origin of system (21) and (23) is based on a modification of the 
mean-field expansion of section 2 with taking into account a particular solution of 
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functional-derivative equation (20). It is easy to see that SDE (20) has the simple 
particular solution 

Z P art(v) = \z2-rf- A -7], (24) 

where functions Z 2 ( X , , ) and A(x — y) satisfy to the system of equations (21) 
V x V J 

and (23). 

This system should be supplemented by one more equation for Z 2 : 

f, l,\zA* y)+xzj* ) z 2 ( [, y 



x y J \ x y J y x y J y x y 

= 5{x y')Z 2 [ ^ ^j+8(x- y")Z 2 ( V yl J . (25) 

The system of three equations (21), (23) and (25) for two functions A and Z 2 
are overfull and seemingly has not physically meaningful solutions. However, as it 
will be show below, equation (25) does not play a role for the construction of the 
modified mean-field expansion. Really, for this construction we define the functional 

G as 

Z = Z part (r ] )+\ogG. (26) 

Taking into account equations (21) and (23) (but not (25)!), we obtain the SDE for 
the functional G: 

. 2 ^oJG 5 2 G 8G SG SG , SG 

m 2 - 2 — - A—- + A— A + r/ ■ — A— Z 2 - r] — XZ 2 ■ r\ — = 

Srj SrjSrj Si] Si] Srj Si] 



= G 



\Z 2 ■ T) Z 2 • T) — T] • Z 2 • T] 



(27) 



(in operator notation; the central dot is an integration) . If we require also to satisfy 
equation (25), then rhs of equation (27) is zero, and we obtain the situation, which 
in the theory of ordinary differential equations is named as " the reduction of order" : 
instead of the second-order linear equation (3), we have a first-order linear equation 
for SG/Sr). But the construction of the mean-field expansion for G, which is similar 
to the mean- field expansion for G of section 2, absolutely has not need of such 
reduction. For this reason we shall not require a satisfaction of equation (25). 

To construct the mean-field expansion for G we consider as a leading-order 
approximation (in full correspondence with the principle of the construction of 
mean- field expansion of section 2) the terms of equation (27), which do not contain 
the source i] manifestly, i.e. the leading-order equation is 

I 2 q2\ . S 2 G . SGq 

{m - d ^H^) - " s^s^x) + x SM^j A{x - y) -° (28) 

All other terms (including the rhs of (27)) are considered as the perturbation. 
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This equation also has a linear exponent as a solution. Then, in correspondence 
with section 2, we have the iteration scheme 



G = G<°> + G« + • • ■ 
and the equation for is 

SG r , , S 2 G n , (5G„ . f SG n - 



(m 2 -d 2 



5rj 5rj5rj 5rj 



1 . <5G„_i <5G„_i 1 

--A — Z2-r)-\Z2-ri — j > 

or/ Or/ J 



Z2 ■ V Z 2 ■ r\ - rj ■ Z2 ■ r\ 



G n 



The general solution of this equation is the functional G^ = P^G^°\ where 
is a polynomial on 77. At nth step of this iteration scheme, we have a closed system 
of linear integral equations, and therefore this scheme is much less complicated in 
comparison to the scheme of preceding section. Equations (21) and (23) are the 
basic approximation for this expansion. 



5 The renormalized equations of the two-particle 
approximation 

Equations (21) and (23) are the system of non-linear equations for the functions 
A and Z2. In equation (23), the two-particle function Z2 can be considered as a 
functional of A, and the "solution" of this equation can be easily found: 

x y 



(29) 



Z 2 ^ x , y yl )=A c (x-y')A(x'-y)- 
- J dxidx 2 A c (x - x 1 )A(x' - x 2 ).f(x 1 - x 2 )A(x 1 - y)A c (x 2 - y'). 
Here, function f(p) in the momentum space is a solution of the equation 

where; 

1 



(30) 



A -(P) = 2j _ 2 - ( 31 ) 

Taking into account equations (29) and (30), we obtain for A the following equation 
in the momentum space 



(m 2 +p 2 )A(p) = l-A(p) J -0^A c (p-q)f(q). 



(32) 
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Note that a very crude approximation for equation (32) is 

(m 2 + P 2 )A(p) » 1 => A(p) » A c (p). (33) 

Within this approximation we obtain as a solution of equation (30) the mean-field 
amplitude of equation (8). So the mean-field approximation and the equivalent 
leading-order l/TV-expansion are contained in the two-particle approximation as 
an approximate solution. 

The renormalization of equations (30) and (32) can be performed in correspon- 
dence with the general recipe of section 2 by introducing counter-term Lagrangian 
(11) and full Lagrangian (12). In terms of the bare quantities the bare propagator 
is 

AjV) = m 2 b +\ b A b (x = 0)+p 2 + X b (p 2 ), (34) 



1 t 

where 



b{p2) = j (^ Ac (P-«)M«) (35) 



is the bare mass operator. 

The normalization conditions (14) and (15) for the renormalized propagator 
A = A b define the mass- renormalization counter-term 8m 2 and the field-renormalization 
constant 

^(l + SUO))- 1 . (36) 
Then the equation for the renormalized propagator is 

(m 2 + p 2 )A(p 2 ) = 1 - A(p 2 )S r (p 2 ), (37) 

where; 

S r (p 2 ) - z[^ b (p 2 ) - E 6 (0) - p 2 K(0)} (38) 

is the renormalized mass operator. 
The equation for f b is 

where 

Lb(p) = / T^yi A ^ + Q) A >>(Q) = z J -J^j A C (P + q)A(q) (40) 

is the bare loop operator. 

The renormalized amplitude F is defined as an amputation of the connected 
part Z% on of the renormalised two-particle function Zi: 

z con = _a • A • F • A • A. 

In correspondence with the bare version of equation (29) 

Z c 2 on = z- 2 Z% n = -z- 2 A c ■ A b ■ f b ■ A b ■ A c = - A c • A • f b ■ A ■ A c 
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and we have (in momentum space): 



P J P J ) = A-\p x )A c (p x )f b (p x + Py )A c (p' y )A-\p' y ). 

fx fy / 

The renormalized coupling A is defined as 

F \ n = A > 
I pi =0 

and due to the normalization condition (14) for the propagator we have 

A(0) = A. (41) 

This normalization condition together with equation (39) defines the coupling- 
renormalization constant z\ and the renormalized equation for f = f b : 

m^J +LAp% (42) 

where 

L r (p 2 ) = L b (p 2 ) - L b (0) (43) 

is the renormalized loop operator. 

Equations (37) and (42) are the system of nonlinear integral equations for the 
propagator and the amplitude. This system with taking into account the normal- 
ization conditions (41) and (14)-(15) can be solved by the expansion in the vicinity 
of the point p = 0. Such solution, however, is not interesting because it is some 
part of the usual perturbation theory over the coupling A. 

Much more interesting problem is to look for the asymptotic behavior of the 
solution at large momenta. In the large-momenta region, an essential technical 
simplification is possible, namely, one can replace in integrals (35) and (40) the 
function A c (see (31)) by massless function 1/p 2 : 

<% 2 ) A c (p - q) =► / — ij — (44) 



(2tt)4 ^ ' ^ *> J (2tt) 4 ip-q) 
Then it is possible to use the well-known formula 

f d 4 q $(q 2 ) 1 [1 [ p \ t 2x 2j 2 [°° ^ 2SJ 



(45) 



The approximation (44) is quite usual in investigations in the deep-Euclidean region, 
and formula (45) highly enables the calculations and, as a major point, permits us 
to go from integral equations to differential ones (see below) . 

As an example of using formula (45), we calculate the field-renormalization 
constant z. Due to the normalization condition (41) we have / — >• A at p — > 0. 
Consequently, with the application of formula (45) we obtain 

*M = ~\-^ (46) 
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and, taking into account (36), 

z= Y^h^- (47) 

Due to the positivity of z, this formula implies the limitation A < 32tt 2 w 320 on 
all following calculations. 

Further essential simplification of the system of equations for / and A consists 
in replacing of equation (37) by the approximate relation 

(m 2 +p 2 )A(p 2 ) « 1 - A c (p 2 )£ r (p 2 ). (48) 

This approximation enables us "to unleash" the system, keeping at the time the non- 
linearity, i.e., to obtain the non-linear integral equation for /, which does not contain 
A. This approximation, of course, is less substantiated and should be considered 
as an iteration of the initial equations. Some grounds for this approximation can 
be obtained a posteriori if the deep-Euclidean behavior of the propagator A does 
not essentially differ from the behavior of the free propagator. 

Approximation (48) defines the equation for /, which will be a main object of 
the following consideration. For the future convenience we introduce new quantities 
y, g and dimensionless variable t as 



m 



The normalization condition for y(t) is 

2/(0) - g- (50) 
After a simple calculation, we obtain the equation for y(t) in the form 



y 

where 



- = -+l(t)+ f dhKit^yih), (51) 
V 9 Jo 



l(t) - (f - l) Iog(l +*) + (!- .9) (1 - \ log(l + *)) , 



(52) 



and the kernel K is 



K(t|tl) = ^_ 1 + I log i±l + (llog ffl±!|. m 

Integral equation (51) is the non- linear Volterra equation and can be reduced to 
the non-linear differential equation of the fourth order 

£[* +1 >"£(;)]-«- J -* < 54 > 

It worth noting that the differential equation (54) will be used only for a suggestion 
the form of the asymptotic solution of the integral equation (51). 
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6 Asymptotic solution for the amplitude 

As it easy to see, the differential equation (54) at g ^ 2 has the exact solution 

Vexact =5-2- (55) 

This solution is not, however, the exact solution of the integral equation (51), 
but the leading term of the asymptotic expansion at t — > oo. Indeed, a simple 
calculation demonstrates that after the substitution of (55) into (51) the leading 
increasing logarithmic term is cancelled. Hence, the asymptotic solution of the 
integral equation (51) at g ^ 2 has the form 

y = g-2 + <p(t), (56) 

where <p(t) — >• at t — >• oo. The function ip(t) at large t satisfies to the linear 
differential equation 



d 2 r ^d 2 



dt 2 



= (.9 - 2)V (57) 



This equation suggests the form of next-to-the- leading asymptotic term ip. The 
solutions of this Euler-type equation are t a_1 , where 

a 2 (a 2 -l) = (g-2) 2 . (58) 

This characteristic equation has two real and two imaginary roots. The first real 
root is greater than 1, and corresponds to an increasing solution. The second real 
root is less than -1, and corresponds to a fast-decreasing solution. The imaginary 
roots correspond to decreasing as t' 1 and oscillating solutions. Hence, just these 
roots define the next-to-the-leading asymptotic term of the solution of integral 
equation (51). 

For further calculation it is convenient to go to the variable 

x = log(l + t) (59) 
In terms of this variable the next-to-the-leading asymptotic term is 

ip = e~ x (Acosujx - Ssinwx), (60) 

where 



VV C9 - 2)2+ i-^ (61) 

Substituting (56) and (60) into the integral equation (51) gives us at t — > oo the 
relation 

- = --l-a(u})A + ub(u})B, (62) 
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where 

1 1 

°H=E (fc + 1)^+^2). fo H = E fc(fc + 1)(fc2+w2) - (63) 

These sums can also be expressed in terms of the special functions (see, e.g. [13]). 

Due to the normalization condition (50) we have A = 2, and equation (62) 
defines the quantity B{g). 



7 Asymptotic behavior of the amplitude 



At t — > oo we have y ^ g—2. At g < 2 such behavior contradicts to continuity. If the 
function y is continuous, then at g < 2 y(to) = in some point to € (0; +oo) (since 
y changes the sign). But in the case the integral equation (51) is not fulfilled in the 
point to- Consequently the function y is singular in some point of pre- asymptotic 
region (Landau pole, or something similar), and at g < 2 this model is inconsistent 
beyond the region of small momenta. 

Hence, the region of applicability of this model is g > 2, or 



where 



A rr <c A <c 2A r 



\ cr = 16tt 2 160. 



(64) 
(65) 



The first two terms of asymptotic expansion in this interval of A at large t = ^ 
are given as follows 



where 



and 



/ ~ 2(A - A cr ) + 

LO(X) 

B(X) = 



2A, r -A/ B(X) 

cos(uj log t) — sin(w log t) 



4(A - A cr 
2A cr A 

(2A cr — A) 2 
4A(A - A cr ) 



At A -> A c 



/ ~ 2(A - A cr ) + 



C 



1 1 

+ 1 + 2a(w) 
1 



(A — A cr ) 2 t 



- sin(wlogi), 



(66) 

(67) 

(68) 
(69) 



where C = (2-7r) 6 /6(0), i.e., in the vicinity of the critical point A cr the amplitude of 
oscillations tends to infinity and the solution is destroyed. 
On the other end of the interval at A — > 2A rr 



/ ~ 2A cr — 4n\/2\ cr — A - sin(wlogt), 



(70) 
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i.e., f(t) —¥ /(O) - the amplitude at t — > oo tends to the initial value of the point 
t = 0. In other words, at A = 2A cr the asymptotic reconstruction of the amplitude 
occurred. 

At A > 2A cr the method of solution is, strictly speaking, cannot be applied due 
to the negative value of the field- renormalization constant z, but if one formally 
assumes such values, the result of equation (66) can be continued to the region of 
strong coupling. 

8 Conclusions 

The most interesting feature of the result obtained is, of course, the existence of the 
critical point A cr , which divides the whole set of coupling values in two region: the 
weak-coupling region of inconsistency and the strong-coupling region of non-trivial 
self-consistent behavior. 4 (The existence of the second specific point seems to be an 
artefact of the calculation scheme.) This division of the coupling values cannot be 
indicated by the expansions tightly connected with perturbation theory (e.g., the 
rcnormalization-group summation and the 1 /iV-expansion) . On the other hand, the 
strong-coupling limit is also not susceptible to the existence of such division, since 
in the framework of the Schwinger-Dyson formalism the strong-coupling expansion 
requires for the physically-acceptable interpretation some subsidiary summation 
[14, 15] based on unknown grounds. 
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